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Abstract 



We present the results of a search for v^, oscillations in the NOMAD experiment 
at CERN. The experiment looked for the appearance of in a predominantly 
wide-band neutrino beam at the CERN SPS. No evidence for oscillations was found. 
The 90% confidence limits obtained are Am^ < 0.4 eV^ for maximal mixing and 
sin^(2^) < 1.4 X 10~^ for large Am^. This result excludes the LSND allowed region 
of oscillation parameters with Am^ > 10 eV^. 

Key words: Neutrino oscillations 
PACS: 14.60.Pq 



1 Introduction 

The NOMAD experiment was designed to search for Ur appearance from neu- 
trino oscillations in the CERN wide-band neutrino beam produced by the 
450 GeV proton synchrotron (SPS). The detector was optimized to identify 
efficiently electrons from — > e^UgUT- decays and therefore could also be 
used to look for appearance in a predominantly beam by detecting their 
charged current (CC) interactions e~X. The main motivation for this 

search was the evidence for z/^ — > z/g and z/^ z/g oscillations found by the 
LSND experiment [1]. For z/^ z/g oscillations with Am^ > 10 eV^ and with 
the probability of 2.6 x 10~^ observed by LSND, a signal should be seen in the 
NOMAD data. The sensitivity of the NOMAD experiment to lower values of 
Am^ is limited by its L/E^ ratio of ~ 0.025 km/GeV, where L is the average 
source to detector distance and E^, is the average neutrino energy. 

Preliminary results of the search for i/^u — > oscillations in NOMAD were 
presented earher [2]. In this paper we report the final results of our "bfind" 
analysis. 



2 NOMAD detector and data collection 

The NOMAD detector [3] consisted of a large dipolc magnet delivering a 
field of 0.4 T and housing several subdetectors, starting with an active tar- 
get composed of 132 planes of drift chambers (DC) of 3x3 m^ [4]. The walls 

* Corresponding author. 

Email address: Slava.Valouev@cern.ch (V. Valuev). 
^ Deceased 
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of the chambers provided a low average density (0.1 g/cm^) target with a 
mass of 2.7 tons. The density of the chambers was low enough to allow an 
accurate measurement of the momenta of the charged particles produced in 
the neutrino interactions. The chambers were followed by nine transition ra- 
diation (TRD) modules [5] each consisting of a polypropylene radiator and 
a plane of straw tubes operated with an 80% xenon and 20% methane gas 
mixture. An electromagnetic calorimeter (ECAL) consisting of 875 lead glass 
blocks [6,7] provided a measurement of the energies of electrons and photons 

with a resolution of a{E)/E = 3.2%/^E{GeV) + 1%. The ECAL was pre- 
ceded by a lead-proportional tube preshower for better photon localization. 
A hadron calorimeter (HCAL) was located just beyond the magnet coil and 
was followed by two muon stations consisting of large area drift chambers, the 
first station located after 8, and the second one after 13 interaction lengths of 
iron. Two planes of scintillator counters, Ti and T2, were placed before and 
after the TRD modules. A third scintillator plane, V, placed upstream of the 
magnet, was used to veto interactions caused by incoming charged particles. 
The trigger [8] used in this analysis was V x Ti x T2. 



3 Neutrino beam 

The CERN West Area Neutrino Facility (WANE) neutrino beam [9] was pro- 
duced by impinging 450 GeV protons extracted from the SPS onto a target 
consisting of beryUium rods adding up to a total thickness of 110 cm. The 
secondary particles emerging from the target were focused into a near parallel 
beam by two magnetic lenses (the horn and the reflector) providing toroidal 
magnetic fields. When running in neutrino mode positively charged particles 
were focused. When running in the antineutrino mode the polarity of the two 
lenses was reversed thus focusing negatively charged particles. The focused 
particles then traversed a 290 m long decay tunnel followed by an iron and 
earth shield. Neutrinos originating from the decay of these particles travelled 
on average a distance of 625 m before reaching the NOMAD detector. 

Since the oscillation search implies a direct comparison between the measured 
and expected ratios of the number of z/g CC to z/^ CC interactions, an accurate 
prediction of the neutrino fluxes and spectra is crucial. They are computed 
with a detailed Monte Carlo simulation of the neutrino beam, referred to as 
NUBEAM and described in Ref. [10]. This is implemented in three steps. First, 
the yields of the secondary particles from the interactions of 450 GeV protons 
with the Be target are calculated with the 2000 version of FLUKA [11], a 
generator of hadronic interactions. These yields are then modified in order to 
agree with all measurements currently available in the relevant energy and 
angular range, namely the SPY/NA56 [12] and NA20 [13] results. Finally, the 
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Fig. 1. Composition of the and v^, energy spectra at NOMAD, within a transverse 
fiducial area of 260 x 260 cm^, as predicted by the NOMAD simulation of the 
neutrino beam line. 

propagation of the secondary particles is described by a simulation program 
based on GEANT3 [14], which includes an accurate description of the magnetic 
field in the horn and reflector, and the modelling of reinteractions in the beam 
elements. 



The resulting energy spectra of i/^ and i^e , and of their components, are shown 
in Fig. 1. The z/^ flux is predominantly due to decays of vr"*" up to 60 GeV 
neutrino energy and to those of beyond this energy. The bulk of the 
flux comes from the decays of K~^, with contributing at the level of about 
18% and fi^ at the level of about 14%. The composition of the beam is shown 
in Table 1. 

Table 1 

Average energies and relative abundances of the fluxes and charged current events 
of the four principal neutrino flavours at NOMAD, within a transverse fiducial area 
of 260 X 260 cm^. 



Flux CC interactions 



Flavour 


{Eu) [GeV] 


Rel. abund. 


(E) [GeV] 


Rel. abund. 




24.3 


1.0 


47.5 


1.0 




17.2 


0.068 


42.0 


0.024 




36.4 


0.010 


58.2 


0.015 


l>e 


27.6 


0.0027 


50.9 


0.0015 



An alternative method of predicting the Uf, content of the beam, not used in 
this paper, has also been studied. This method, referred to as the Empirical 
Parameterization, does not use hadronic interaction packages such as FLUKA 
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to predict the yield of particles from p-Be interactions. Instead it uses the 
CC, z/^ CC and z/g CC events observed in NOMAD as well as the results of 
NA20 and NA56 to constrain the hadron production cross sections, and using 
these, to predict the content of the beam. 



4 Event simulation 

The neutrino flux generated by NUBEAM was used as an input to the NO- 
MAD event generator to produce interactions of u^, u^, v^. and v^- Deep- 
inelastic scattering (DIS) events were simulated with a modified version of 
the LEPTO 6.1 event generator [15], with and W'^ cutoffs removed. Quasi- 
elastic (QEL) [16] and resonance production (RES) [17] events were generated 
as well. The GRV-HO parameterization [18] of the parton density functions 
and the nucleon Fermi motion distribution of Ref. [19], truncated at 1 GeV/c, 
were used along with JETSET 7.4 [20] to treat the fragmentation. 

The secondary particles produced in these interactions were then propagated 
through a full GEANT3 [14] simulation of the detector. The size of the simu- 
lated samples exceeded the data samples by about a factor of three for CC, 
10 for CC and neutral current (NC) and 100 for v^, CC interactions. 

The contributions of QEL, RES and DIS events in the Monte Carlo were 
adjusted to reproduce the distribution of CC interactions observed 
in the data, taking into account the non-isoscalarity of the NOMAD target 
(52.4% protons and 47.6% neutrons). The fractions of QEL and RES events 
included in the Monte Carlo are hsted in Table 2. 

Table 2 

Estimated percentages of QEL and RES events in the NOMAD data. 



Vf. f e 



QEL 


RES 


QEL 


RES 


QEL 


RES 


QEL RES 


2.5% 


3.3% 


6.3% 


7.1% 


1.7% 


2.2% 


4.3% 4.8% 



5 Data collection and analysis 

NOMAD collected data from 1995 to 1998. Most of the running, a total expo- 
sure of 5.1 X 10^^ protons on target (pot), was in neutrino mode and yielded 
1.3 X 10^ CC interactions in the fiducial volume of the detector. Some 
data, amounting to 0.44 x 10^^ pot, were also collected in antineutrino mode 
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and some, 0.04 x 10^^ pot, in zero-focusing mode (with the horn and reflector 
switched off); these data were used mostly to check the beam simulation. 

The trajectories of charged particles are reconstructed from the hits in the 
drift chambers and, from these trajectories, momenta are computed using 
the Kalman filter technique [21], which accounts for energy loss along the 
trajectory. As a first step the energy loss model used is that for pions, resulting 
in a momentum estimate, p^^, at the beginning of the track. Particles later 
identified as electrons or positrons arc refitted [22] with an additional average 
energy loss due to bremsstrahlung, resulting in a new estimate, Pe, of the 
momentum. Energy clusters in the ECAL not associated to charged particles 
are assumed to be due to photons. 

Vertices are reconstructed from the trajectories of charged particles. The en- 
ergy of the incident neutrino, E„, is approximated by the total (visible) energy 
of an event computed from the sum of the energies of all observed primary 
particles and of photons. The reconstruction efficiency for the hadronic jet 
was found to be overestimated by the Monte Carlo. This was due mostly to 
reconstruction effects such as the density cut described in Section 7, as well 
as an inadequate treatment of nuclear reintcractions (the interactions of pro- 
duced particles inside the nucleus in which the neutrino interaction occurred) 
and a harder fragmentation in the Monte Carlo than observed in the data. 
The resulting overestimate of the scale of the hadronic energy. Eh, could be 
studied by noting that the differential cross section for charged current neu- 
trino interactions is almost independent of Ibj = EfjEi,. This entails that the 
distribution of Re, the ratio of the average neutrino energy in the data to the 
average neutrino energy in the Monte Carlo, will be independent of Ibj if Eh is 
correctly measured. However, if the measured E^, is systematically reduced by 
an amount a with respect to the true value of E^, then the dependence of Re 
on Ybj can be described by a straight line with slope (1 — a) for small values 
of a. Fitting the Re distribution as a function of Ibj allowed us to determine 
that the correction to Eh is a — (8.3 ± 1.5)%. 

Since the electron radiates bremsstrahlung photons in traversing the drift 
chambers, in order to have an accurate measure of its energy, these photons 
must be identified and their energy added to the energy of the ECAL cluster 
at the end of the electron trajectory. Because of the curvature of the electron 
trajectory in the magnetic field these photons are located in a vertical fan 
delimited, on the one hand, by the actual trajectory of the electron between 
the event vertex and the point of impact of the electron on the ECAL and, on 
the other hand, by the extrapolation of the initial direction of the electron to 
the ECAL. The energy of photons in the ECAL and of photon conversions in 
the DC found in this region is included, resulting in a measure of the electron 
energy, E'brem, with an average resolution of 2.1%. 
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6 Principles of oscillation search 

The Vfj, — > Ue oscillation signal should manifest itself as an excess in the number 
of Vf, CC events over that expected for an intrinsic contamination in the 
beam (about 1% of z/^). In order to reduce systematic uncertainties associated 
with absolute flux predictions and selection efficiencies, we study the ratio 
Ren of the number of to charged current interactions. Due to different 
energy and radial distributions of incident electron and muon neutrinos, the 
contribution of the intrinsic i/g component is smaller at low energies, E^, 
where a low Am^ oscillation signal is expected, and at small radial distances 
from the beam axis, r. Thus, the sensitivity of the search is increased by taking 
into account the dependence of Re^ on Ei, and r. 

The presence or absence of — > Ug oscillations is established by comparing 
the measured i?e/i with the one expected in the absence of oscillations. In 
order to avoid biases, we adopted a "blind analysis" strategy: the comparison 
of the measured to the predicted R^^ is not made until the accuracy of the flux 
predictions and the robustness of the data analysis have been demonstrated 
and until all selection criteria are flxed. A number of control data samples in 
which no oscillation contribution is expected (charged current interactions of 
Ufj,, and z/e in neutrino mode, and of z/^, i/^, Ug and Ug in antineutrino and 
zero-focusing modes) are used to verify the flux predictions and the validity of 
the Monte Carlo simulation [10]. It should be noted that no oscillation signal 
is expected to be measurable in z/g since the intrinsic ratio of Uf_i/ue in the 
beam is four times smaller than the intrinsic V/j,/ Ve ratio and the antineutrino 
statistics is limited. 



7 Event selection 

In order to calculate R^.^, pure samples of z/g CC and z/^ CC interactions are 
selected. The initial data sample for i/g CC interactions is complementary to 
that used in the Un CC selection described below, i.e., it consists only of those 
events that include no muon (identifled with looser criteria than in the CC 
selection). The basic requirement is the presence of a track associated with 
the neutrino interaction vertex, pointing to an energy deposition in the ECAL 
and identified as an electron in the TRD and ECAL. The identification criteria 
are: 

• Pulse heights in the TRD consistent with those of an electron and such that 
isolated charged pions are rejected by at least a factor of 1000. 

• A momentum-energy match satisfying: 

• (-Ebrem -P7r)/(£^brem+P7r) > "O-S; 
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• (£^brem - Pe)/(£^brcm + Pe) < 0.4; 

• no activity in HCAL associated with the electron trajectory. 
Electrons from conversions and Dalitz decays are rejected by requiring: 

• that the first point on the candidate track be within 15 cm of the primary 

vertex; 

• that no positively charged track, either identified as a positron in the TRD 
or missing the TRD, when taken together with the candidate electron, re- 
sults in the combination being consistent with a conversion. The criteria 
used are the invariant mass and the angle between the plane containing the 
trajectories of the two tracks and the vertical. 

In order to reduce further the background from neutral current and charged 
current events in which the muon was not identified, kinematic cuts are also 
applied using the following two variables: 

• (peh, the angle between the electron and the hadronic jet in the plane trans- 
verse to the neutrino beam direction. 

• Qiep, the component of the electron momentum perpendicular to the had- 
ronic jet direction. 

Charged pions simulating secondary electrons and conversion electrons are 
part of the hadronic jet, resulting in small values of these variables. Primary 
electrons are isolated resulting in large values of 0eh and Qiep- These differences 
are evident in the 0eh-Qiep plots shown in Fig. 2 on which the cut used is also 
shown. 

Single track events, originating mostly from quasi-elastic and resonance pro- 
duction interactions, are treated somewhat differently because they are more 
prone to background arising from charged particles entering the detector with- 
out being registered in the veto counter and because no hadronic jet containing 
primary charged particles is present. The criteria used to select them are: 

• no activity in the drift chambers upstream of the beginning of the track and 
consistent with having given rise to the track; 

• the angle between the single electron and neutrino beam directions, 9i,e, 
smaller than 150 mrad; 

• l^Bj < 0.5; 

• PeX > 2me, with mg the mass of the electron, a cut that rejects neutrino- 
electron scattering events. 

Lastly, only events with pe > 2.5 GeV/c and < 300 GeV are retained. 
These selection criteria result in an efficiency for i/g CC, estimated from the 
Monte Carlo, of 43.9%. 
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Fig. 2. The two-dimensional distributions (/!>eh-Qiep (defined in the text) for Monte 
Carlo CC and NC (background) events and for Vf. CC (signal) events, as well 
as for the NOMAD data. The events to the right of the curve were selected. 



The surviving background contribution to the z/g CC sample is estimated, 
from the Monte Carlo, to be 1.8%. It consists mostly of electrons from photon 
conversions. Their rejection depends critically on the reconstruction of the 
accompanying positron and on identifying a conversion point distinct from 
the primary vertex. The reconstruction of very low energy positrons and the 
separation of a conversion vertex from the primary vertex in a high multiplicity 
event could be different in the data and in the Monte Carlo. As a cross-check, 
a class of e+ events that failed the kinematic cuts or were produced far from 
the primary vertex, and thus consisting almost entirely of background, were 
selected in both the data and in the Monte Carlo. The number of such events 
was found to be higher in the data by a factor that varied between 1.0 and 
1.6 depending on Ey. The Monte Carlo background estimate for the events 
was therefore multiplied by this same factor, thus raising the total background 
estimate to 2.3%. 
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Fig. 3. The distributions of z (defined in the text) for CC (left) and fg CC (right) 
candidates in the data (points with error bars) and in the Monte Carlo (histogram). 
The Monte Carlo distribution of CC events is normalized to the number of CC 
events in the data; the Monte Carlo distribution of Vf, CC events is normalized using 
the relative CC/v^ CC abundance predicted by our simulation. The background 
contribution to the CC sample is shown in the hatched histogram. 

Charged current interactions of are characterized by the presence of a pri- 
mary muon in the final state, which had to penetrate 13 interaction lengths 
of absorber material to reach both muon stations in order to be identified. 
In addition, in order to minimize the differences between selection efficiencies 
of CC and z/g CC events, we apply kinematic criteria identical to those 
used in the Ve CC selection, although they are not needed for the background 
suppression. The resulting CC data sample has a negligible background 
contamination; the average selection efficiency is 60%. 

The geometrical and kinematical distributions of both types of events are 
well reproduced by the Monte Carlo simulation, with the exception of the 
distribution of interaction vertices along z, the beam direction. The data and 
Monte Carlo z distributions are shown in Fig. 3 for CC and Ve CC events. 
It can be seen that fewer events are present in the data than in the Monte 
Carlo at small z especially for CC events. The origin of this difference is due 
mostly to a cut introduced during the reconstruction of events: events with 
a very high density of hits in the drift chambers were not reconstructed due 
to a prohibitive reconstruction time. Since the data has on average a higher 
density of hits than the Monte Carlo, the effect of this cut is different on 
the two samples. Furthermore, since electrons radiate photons in traversing 
the drift chambers and some of these photons convert, the density of hits 
in z/g CC events is large thus exacerbating the effect of the density cut for 
these events. The reprocessing of a sample of data and Monte Carlo events 
without this density cut resulted in z distributions that were in much better 
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Visible energy (GeV) Visible energy (GeV) 

Fig. 4. Neutrino energy spectra for the data (points with error bars) and the Monte 
Carlo (histogram), for i/^ CC (left) and Ve CC (right) candidates. The normaliza- 
tion of the Monte Carlo distributions is the same as in Fig. 3. The background 
contribution to the z^e CC sample is shown in the hatched histogram. 

agreement. We therefore decided to restrict the analysis to events occurring 
in the 72 downstream planes of drift chambers by requiring z > 184 cm. This 
restriction resulted in a loss of about 30% of the Ue CC events and 35% of the 
CC events. It should be noted that any oscillation effects could not manifest 
themselves over this distance since the point of origin of the neutrinos is spread 
over more than 300 m. 

A total of 5,584 CC and 472,378 z/^ CC events were retained. Their energy 
spectra and radial distributions are shown in Fig. 4 and Fig. 5, respectively. 



8 Systematic uncertainties 

The single largest uncertainty in this oscillation search is the uncertainty in 
the prediction of the fraction and the energy spectrum of intrinsic z/g present 
in the beam. The computation of this uncertainty is described in detail in 
Ref. [10]. Most of it is due to the limited knowledge of the particle production 
cross sections in p-Be interactions. In turn, this is due to the small number 
of experimental data points on tt"*" and production measured by NA20 
and NA56, especially at non-zero values of transverse momentum. This uncer- 
tainty is energy- dependent; its typical fractional value at low is 4%. The 
second largest systematic error of about 3% comes from an uncertainty in 
the production of which accounts for 18% of the flux. Other potential 
sources of errors (such as tertiary particle yields, variations in the horn current, 
misalignments of the focusing devices and collimators, or inaccuracies in the 
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Fig. 5. The distribution of r^, the square of the radial position of the neutrino 
interaction vertex with respect to the nominal beam axis, for the data (points with 
error bars) and the Monte Carlo (histogram), for CC (left) and v^, CC (right) 
candidates. The normalization of the Monte Carlo distributions is the same as in 
the two previous plots. The background contribution to the CC sample is shown 
in the hatched histogram. 

simulation of the beam line elements) have also been investigated [10]; their 
cumulative contribution is about 3%. The overall uncertainty arising from the 
knowledge of the beam composition is divided into an energy-independent, or 
normalization, uncertainty and an energy-dependent one. The normalization 
uncertainty on R^.^ is 4.2%, while the energy- dependent uncertainty, shown in 
Fig. 6, varies from 4% to 7%. 
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Fig. 6. Energy-dependent uncertainty in the prediction of the Ren ratio. 
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The following additional sources of systematic uncertainties arising from the 
data selection and analysis were studied: 

• The correction to the background contribution. By studying the e"*" events 
used to determine this correction, the uncertainty attached to it was found 
to be 15%, resulting in an uncertainty of less than 1% in 

• The electron reconstruction and identification. The electron reconstruction 
efficiency in the drift chambers (~ 98%) was studied by defining electron 
tracks using the TRD, preshower and ECAL only, and then computing the 
frequency that a DC track is associated to the electron track. The elec- 
tron identification efficiency of the TRD (~ 93%) was studied using 5-rays 
produced by isolated muons arising from a nearby test beam and travers- 
ing the NOMAD detector outside of the neutrino spill. The efficiency of 
the momentum-energy consistency check was studied using from photon 
conversions. Taken together, the uncertainties in the reconstruction and 
identification efficiencies for electrons result in a 1% uncertainty in i?e^. 

• A possible relative shift between the electron and muon energy scales. The 
upper hmit on the absolute energy uncertainty of the ECAL was 0.5% [7], 
resulting in an uncertainty of less than 1% in i?e/i- Non-linearity correc- 
tions to the ECAL energy scale have negligible effects on Ren- The muon 
momentum scale introduced a negligible error. 

• The mixture of QEL, RES and DIS events in the Monte Carlo. Cross sections 
for QEL and RES processes were varied by as much as 50% resulting in a 
neghgible effect on Re^. 

• The fragmentation model. The effect on i^e^t of the fragmentation and nu- 
clear reinteraction model used in the Monte Carlo was studied by comparing 
the prediction for i?e/i of two such models. Differences of up to 2% were ob- 
served and introduced as an energy-dependent systematic uncertainty. 

• The kinematic (isolation) cuts. This was studied by tightening and loosening 
the cuts in both the data and Monte Carlo and studying the effect on Re^. 
Since this was a blind analysis, this study could only be made after the 
analysis was frozen and the electron neutrino data examined. The effect 
was negligible in the relevant neutrino energy range, from 10 to 80 GeV. 

• The hadronic energy scale. The 18% uncertainty on a, the hadronic energy 
correction, resulted in less than a 1% uncertainty on Refj, throughout the 
neutrino energy range. It was included as an energy-dependent uncertainty. 
The uncertainty arising from the method of calculating the hadronic en- 
ergy was studied using two alternative methods to calculate it - the Myatt 
method [23] and the double-angle method [24] , which use both the energy 
and direction of the lepton and only the direction of the hadronic jet. Af- 
ter computing the appropriate energy correction for each of these methods, 
no systematic differences were observed between these alternative methods 
and the method described in Section 5. Again, this study could only be 
performed after examining the electron data. 
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9 Results 



The Refj, distribution as a function of the visible energy obtained from the data 
is shown in Fig. 7, for the full radial acceptance (left) and in three radial bins 
(right). It is in good agreement with the Monte Carlo prediction under the 
no-oscillation hypothesis, also shown in the figure as ila uncertainty bands: 
a of 37.1/30 d.o.f. is obtained when the data are analyzed and compared 
to the simulation in the 10 energy bins and the 3 radial bins shown in Fig. 7 
(incorporating both statistical and systematic uncertainties). The best fit to 
^'u oscillations, in the two-family approximation, gives a similar chi- 



squared value. 
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Fig. 7. The Re^ ratio as a function of the visible energy for the data (points) and 
for the Monte Carlo prediction assuming no oscillations (filled bands), for the full 
radial acceptance (left) and in the three radial bins (right). The upper and lower 
boundaries of the bands correspond to the predictions with ztlcr uncertainty, where 
a includes both the normalization and energy-dependent systematic uncertainties 
added in quadrature. 

We use a frequentist approach [25] to set a 90% confidence upper limit on the 
oscillation parameters. The resulting exclusion region is shown in Fig. 8, to- 
gether with results of other accelerator experiments, LSND [1], KARMEN [26], 
CCFR [27] and NuTeV [28], and the combined limit of Bugey [29] and Chooz 
[30] reactor experiments. Values of Am? > 0.4 eV^ for maximal mixing and 
sin^(2^) > 1.4 X 10~^ for large Am^ are excluded. For comparison, the sen- 
sitivity [25] of the experiment is found to be Am^ > 0.4 eV^ for maximal 
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Fig. 8. The 90% C.L. exclusion region in the Am^ - sin^(2^) plane from this analysis 
superimposed on the results of other experiments. 

mixing and sin^(2^) > 1.3 x 10~^ at large Am^. Our result rules out the in- 
terpretation of the LSND measurements in terms of z/^ — > z/g oscillations with 
Am2 > 10 eV2. 

This result is less stringent than our preliminary result [2] because of a bet- 
ter understanding of the systematic uncertainties arising from the knowledge 
of the beam composition. In particular, our previous analysis had underesti- 
mated the uncertainty arising from the contribution to the v^. spectrum. In 
addition a different split between normalization and energy-dependent errors 
was implemented for the uncertainty arising from mesons produced in sec- 
ondary interactions. After opening the box defined for the purpose of the blind 
analysis, no modifications were made either to the central value of the beam 
prediction or, other than the z > 184 cm cut described above, to the event 
selection procedure. We have ensured that including the events at 2; < 184 cm 
in the analysis did not alter our limit significantly. 
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10 Conclusion 



The results of a search for — > neutrino oscillations in the NOMAD 
experiment at CERN have been presented. The experiment looked for the ap- 
pearance of Uc in a predominantly z/^ wide-band neutrino beam at the CERN 
SPS. No evidence for oscillations was found. The 90% confidence limits ob- 
tained are Am^ < 0.4 eV^ for maximal mixing and sin^(2^) < 1.4 x 10~^ for 
large Am^. This result excludes the high Am^ region of oscillation parameters 
favoured by the LSND experiment. 
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